(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property 
Organization 
Internationa] Bureau 

(43) International Publication Date 
21 October 2004 (21.10.2004) 




PCT 



(10) International Publication Number 

WO 2004/089217 Al 



(51) International Patent Classification 7 : 
(21) International Application Number: 



A61B 6/03 



(22) International Filing Date: 

(25) Filing Language: 

(26) Publication Language: 



PCT/IB2004/001061 
2 April 2004 (02.04.2004) 
English 
English 



(30) Priority Data: 
03100975.6 



10 April 2003 (10.04.2003) EP 



(71) Applicant (for DE only): PHILIPS INTELLECTUAL 
PROPERTY & STANDARDS GMBH [DE/DE]; Stein- 
damm 94, 20099 Hamburg (DE). 

(71) Applicant (for all designated States except DE, US): 
KONINKLUKE PHILIPS ELECTRONICS N. V. 
[NL/NL]; Groenewoudseweg 1, 5621 BA Eindhoven (NL). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): KOHLER, Thomas 
[DE/DE]; c/o Philips Intellectual Property & Standards 



GmbH, Weisshausstr. 2, 52066 Aachen (DE). GRASS, 
Michael [DE/DE]; c/o Philips Intellectual Property & 
Standards GmbH, Weisshausstr. 2, 52066 Aachen (DE). 
PROKSA, Roland [DE/DE]; c/o Philips Intellectual 
Property & Standards GmbH, Weisshausstr. 2, 52066 
Aachen (DE). 

(74) Agent: MEYER, Michael; Philips Intellectual Property & 
Standards GmbH, Weisshausstr. 2, 52066 Aachen (DE). 

(81) Designated States (unless otherwise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 
AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, TT, 
GB, GD, GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, KB, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, NI, NO, NZ, OM, PG, 
PH, PL, PT, RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, 
TN, TR, TT, TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM, 
ZW. 

(84) Designated States (unless otherwise indicated, for every 
kind of regional protection available): ARIPO (BW, GH, 

[Continued on next page J 



(54) Title: COMPUTER TOMOGRAPHY METHOD FOR A PERIODICALLY MOVING OBJECT 







































l-H 


S. 2 


< 


S-1 






1— 1 








o\ 




00 
























O 










(57) Abstract: The invention relates to a 
computer tomography method in which a 
periodically moving object, in particular 
an organ of the body, is irradiated by a 
cone-shaped beam cluster (4) along a 
trajectory which runs on a cylindrical 
surface. The radiation transmitted through 
the object is measured by means of a 
detector unit (16), and at the same time 
the periodic movement of the object is 
recorded. In order to reconstruct the 
absorption distribution of the object, the 
measured values or the corresponding 
beams are rebinned to form a number of 
parallel projections, where for each of these 
projections a measured value is determined 
whose beam irradiates the object. The 
point in time at which this measured value 
was acquired is allocated to the respective 
projection. For the reconstruction, which 
may for example be carried out using a 
filtered back-projection, only projections 
whose allocated points in time lie within a 
predefined, specific time range (HI) within 
a period of the object movement are used. 
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Computer tomography method for a periodically moving object 



The invention relates to a computer tomography method in which a 
periodically moving object, in particular an organ of the body, which is located in an 
examination area is irradiated by a cone-shaped beam cluster. The invention also relates to a 
computer tomography scanner for carrying out the method and to a computer program for 
5 controlling the computer tomography scanner. 

Within the context of this invention, the expression "periodic movement" 
refers to movements in which a series of object states are repeatedly assumed by the object, 
always in the same sequence, during a measurement. An object state may in this case be 
defined, for example, by the position of the object in the examination area and the shape of 

10 the object or by a specific measured signal which is allocated to the object state. This 
measured signal may, for example, be the signal of an electrocardiograph if the moving 
object is a heart. If an electrocardiograph supplies the same measured signals at different 
points in time, then it is assumed that the object states at these points in time are the same. 
Within the context of the invention, however, the term "periodic movements" also means 

15 movements which could be referred to as "quasi-periodic" and in which object states that are 
repeatedly assumed by the object are not exactly identical but are substantially identical. The 
same applies to the time duration between two substantially identical, repeating object states, 
which time duration need only be substantially constant during the measurement in order to 
be able to refer to the corresponding movement as "periodic" within the context of this 

20 invention. Different period time durations are "substantially identical" if their time difference 
is small compared to the periodic time durations for the respective application. Accordingly, 
two object states are substantially identical if their difference is small compared to the 
differences which the object passes through within an entire period. Therefore, for example, 
movements of pulsating technical or biological objects such as organs of the body or veins 

25 may be referred to as "quasi-periodic" and hence also as "periodic". Moreover, the term 
"periodic movement" also encompasses attenuated oscillations or movements in which the 
object repeatedly assumes only some, that is to say not all, states at substantially identical 
time intervals. 
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In known methods of the type mentioned above, the spatial profile of the 
absorption or of the attenuation of the radiation in the periodically moving object is 
reconstructed from measured values acquired using a detector unit. In this case, the 
movement of the object leads to the measured values containing information from different 

5 object states, and this leads to movement artefacts in the reconstructed data record. 

Therefore, in known reconstruction methods use is only made of measured 
values which have been acquired in time ranges in which the object has moved relatively 
little, that is to say in which the object states change only slightly. In one of these methods, 
the beams firstly undergo parallel rebinning in order to reduce the computational complexity. 

10 In parallel rebinning, the cone-shaped beam cluster for each radiation source position is split 
on the trajectory into beam fans, the beams of which in each case lie in a plane which runs 
parallel to the central axis of the circular or helix-like trajectory. These beam fans are then 
divided into groups such that each group contains only beam fans that are parallel to one 
another. The beam fans of a group are referred to as a projection. A point in time is allocated 

1 5 to each projection, said point in time being the same as the point in time at which the 

measured values corresponding to the central beam fan of the projection were acquired. It is 
assumed that all measured values of a projection were acquired at this point in time. 

In the reconstruction of the absorption distribution from the rebinned beams or 
measured values, use is only made of projections which were acquired during predefinable 

20 time ranges in which the movement of the object is relatively small. Since - apart from in 
respect of the central beam fan - the point in time of a projection does not coincide with the 
actual acquisition time of a measured value of this projection, for the reconstruction use is 
frequently made of measured values whose beams do indeed come into contact with the 
object but were acquired at points in time which lie outside the predefined time ranges. This 

25 leads to movement artefacts which increase as the distance of the object from the central 
beam fan increases. 

It is an object of the present invention to specify a method in which these 
movement artefacts are less pronounced. This object is achieved according to the invention 
by a computer tomography method comprising the steps 

30 - using a radiation source (S) to generate a cone-shaped beam cluster (4) which 

passes through an examination area (13) and a periodically moving object which is located in 
the examination area (13), 

producing a relative movement between the radiation source (S) on the one 
hand and the object located in the examination area (13) on the other hand, where a 
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trajectory, along which the radiation source moves relative to the object, runs on an 
imaginary cylindrical surface that surrounds the object, 

using a detector unit (16) to acquire measured values which depend on the 
intensity in the beam cluster (4) on the other side of the object, during the relative movement, 
5 - recording the periodic movement of the object during the acquisition, 

reconstructing a spatial distribution of the absorption of the periodically 
moving object from the measured values with the aid of the recorded periodic movement of 
the object, comprising the steps: 

a) determining the spatial area taken up by the object in the examination area 
10 (13), 

b) subjecting the measured values to parallel rebinning in order to form a number 
of groups, where the beams corresponding to the measured values of each group form beam 
fans (41 . . .45) which lie in planes that are parallel to one another and to the axis of rotation, 

c) determining for each group a measured value whose beam irradiates the spatial 
1 5 area taken up by the object, and allocating to the respective group the point in time at which 

this measured value was acquired, 

d) determining those groups whose points in time, allocated in step c), lie within 
periodic, predefined time ranges (Hi), 

e) reconstructing the absorption distribution in the object from the measured 
20 values belonging to the groups determined in step d). 

In contrast to known methods, in the invention there is determined for each 
group or projection a measured value whose beam irradiates the object, and the point in time 
at which this measured value was acquired is allocated to the projection. This leads to the 
time intervals between the point in time allocated to a projection and the acquisition times of 
25 said measured values of this projection, which irradiate the object, being smaller than in 
known methods. This reduces the incorrect use of measured values which lie outside the 
predefined time ranges. Since the time ranges are generally selected such that they correlate 
with phases of the object movement which do not contain much movement, movement 
artefacts are reduced. 

30 Claim 2 describes one refinement in which the spatial area taken up by the 

object in the examination area is determined with a low computational complexity by 
reconstruction and segmentation of the object. 

In Claim 3, for each group there is determined the measured value whose 
beam runs through the geometric center of the object. The point in time at which this 
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measured value was acquired is allocated to the respective group. This leads to a further 
reduction in the incorrect use of measured values which lie outside the predefinable time 
range, and this further reduces the movement artefacts. 

Claim 4 shows one refinement in which the periodically moving object is a 
5 heart and in which reconstructed images of a high quality can be produced inter alia with the 
aid of an electrocardiograph. 

In claim 5, the predefinable time ranges are selected such that the object 
moves less in these time ranges than in other time ranges, and this leads to a further reduction 
of the movement artefacts in the reconstructed images. 
10 Claim 6 defines one refinement in which the reconstruction is carried out 

using a filtered back-projection, and this leads to a good image quality of the reconstructed 
object while entailing a low computational complexity. 

In claim 7 the radiation source moves relative to the examination area on a 
helix-like or circular trajectory, and this leads to reconstruction results with a good image 
15 quality. 

A computer tomography scanner for carrying out the method is described in 

claim 8. 

Claim 9 defines a computer program for controlling a computer tomography 
scanner as claimed in claim 8. 



The invention will be further described with reference to examples of 
embodiments shown in the drawings to which, however, the invention is not restricted. 

Fig. 1 shows a computer tomography scanner which can be used to carry out 
25 the method according to the invention. 

Fig. 2 shows a flowchart of the method according to the invention. 

Fig. 3 shows a beam cluster of parallel rebinned beam fans. 

Fig. 4 shows a flowchart for a back-projection. 

Fig. 5 shows an electrocardiogram. 
30 Fig. 6 shows a flowchart for a further back-projection. 

Fig. 7 shows the dependence of an a factor on time. 

Fig. 8 shows the dependence of a p factor on time. 
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The computer tomography scanner shown in figure 1 comprises a gantry 1 
which can rotate about an axis of rotation 14 that runs parallel to the z direction of the 
coordinate system shown in figure 1 . For this purpose, the gantry 1 is driven by a motor 2 at a 
preferably constant but adjustable angular velocity. A radiation source S, for example an 
X-ray generator, is attached to the gantry 1. Said radiation source is provided with a 
collimator arrangement 3 which masks, from the radiation generated by the radiation source 
S, a cone-shaped beam cluster 4, that is to say a beam cluster which has, both in the z 
direction and in a direction perpendicular thereto (i.e. in a plane perpendicular to the axis of 
rotation), a finite dimension other than zero. 

The beam cluster 4 penetrates a cylindrical examination area 13 in which there 
is located a periodically moving object (not shown). In this example of embodiment, said 
object is a beating heart which carries out intrinsic movements and in some circumstances 
also moves back and forth on account of the patient's breathing movements. In other 
embodiments, other periodically moving organs of the body, such as the liver or the brain, 
periodically moving parts of organs of the body or periodically moving technical objects 
could also be irradiated. 

After passing through the examination area 13, the beam cluster 4 comes into 
contact with a detector unit 16 which is attached to the gantry 1 and has a detector surface 
comprising a large number of detector elements which in this embodiment are arranged in 
rows and columns in the form of a matrix. The detector columns preferably run parallel to the 
axis of rotation 14. The detector rows are located in planes perpendicular to the axis of 
rotation, in this embodiment on an arc of a circle around the radiation source S (focus- 
centered detector). However, in other embodiments they may also be formed differently, for 
example describe an arc of a circle around the axis of rotation 14 or be rectilinear. Each of 
the detector elements with which the beam cluster 4 comes into contact supplies, in each 
position of the radiation source, a measured value for a beam from the beam cluster 4. 

The spread angle of the beam cluster 4, designated a max , determines the 
diameter of the object cylinder within which the object to be examined is located during 
acquisition of the measured values. The spread angle is defined as the angle enclosed by a 
beam lying in a plane perpendicular to the axis of rotation 14 at the edge of the beam cluster 
4 with a plane defined by the radiation source S and the axis of rotation 14. The examination 
area 13 or the object or the patient table may be moved parallel to the axis of rotation 14 or to 
the z axis by means of a motor 5. However, the gantry may also be moved in this direction in 
an equivalent manner. When a technical object is concerned rather than a patient, the object 
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can be rotated during an examination while the radiation source S and the detector unit 16 
remain stationary. 

With the aid of the motors 2 and 5, the radiation source S and the detector unit 
1 6 can describe a trajectory relative to the examination area 13, said trajectory running on an 
imaginary cylindrical surface. This trajectory may, for example, run in a helix-like manner 
when both motors are operating. If, on the other hand, the motor 5 for advance in the 
direction of the axis of rotation 14 is idle and the motor 2 is making the gantry rotate, the 
result is a circular trajectory for the radiation source S and the detector unit 16 relative to the 
examination area 13. In this example of embodiment, only the helix-like trajectory is 
considered. 

During the acquisition of the measured values, the movement of the heart is 
recorded in a known manner by means of an electrocardiograph 8. For this purpose, the chest 
area of a patient is connected to the electrocardiograph 8 by way of electrodes (not shown). 
In other embodiments, in particular in the case of other moving objects, the movement of the 
object may be traced in other ways. Thus, for example, the movement information could be 
obtained from the values measured by the detector unit 16 themselves, so that there is no 
need for the movement to be recorded using an additional device, such as an 
electrocardiograph. For this, firstly a kymogram is produced from the measured values, from 
which kymogram the movement can be derived in a known manner. A detailed description of 
this method can be found in "Kymogram detection and kymogram-correlated image 
reconstruction from subsecond spiral computed tomography scans of the heart," M. 
KachelrieG, D.A. Sennst, W. Maxlmoser, W.A. Kalender, Medical Physics 29(7): 1489-1503, 
2002, to which reference is hereby made. 

In this example of embodiment, it is assumed that the patient is not breathing 
during the measurement. The breathing movement can thus be disregarded. In other 
embodiments, the breathing movement could be measured, for example, using a deformable 
chest strap which is connected to a breathing movement measurement device. 

The measured values acquired by the detector unit 16 are fed to a 
reconstruction and image processing computer 10 which is connected to the detector unit 16,. 
for example via contactless data communication (not shown). Moreover, the 
electrocardiogram is transferred from the electrocardiograph 8 to the reconstruction and 
image processing computer 10. The reconstruction and image processing computer 10 
reconstructs the absorption distribution in the examination area 13 and displays it on a 
monitor 1 1 for example. The two motors 2 and 5, the reconstruction and image processing 



WO 2004/089217 PCT/IB2004/001061 

7 

computer 10, the radiation source S, the electrocardiograph 8 and the transfer of the measured 
values from the detector unit 16 to the reconstruction and image processing computer 10 are 
controlled by the control unit 7. The control unit 7 furthermore controls the transfer of the 
electrocardiogram from the electrocardiograph 8 to the reconstruction and image processing 
5 computer 10. 

In other embodiments, the acquired measured values and the measured 
electrocardiograms can be fed for reconstruction firstly to one or more reconstruction 
computers which forward the reconstructed data e.g. via a fiberoptic cable to the image 
processing computer. 

10 Figure 2 shows the sequence of a measurement and reconstruction method 

which can be carried out using the computer tomography scanner shown in figure 1 . 

After initialization in step 101, the gantry rotates at an angular velocity which 
in this example of embodiment is constant. However, it may also vary, for example as a " 
function of time or of the radiation source position. 

15 In step 103, the examination area or the object or the patient table is moved 

parallel to the axis of rotation and the radiation of the radiation source S is switched on so 
that the detector unit 16 can detect the radiation from a large number of angular positions. At 
the same time as or before the switching on of the radiation source S, the electrocardiograph 
8 is activated so that an electrocardiogram is measured at the same time. 

20 Thereafter, the spatial area in which the heart is located within the examination - 

area is determined. For this, the absorption distribution in the examination area is firstly 
reconstructed from the measured values at a low resolution without taking the 
electrocardiogram into account. There is a low resolution for example when a volume of 20 x 
20 x 20 cm 3 is represented by 64 3 voxels. 

25 In step 105, the measured values are parallel rebinned for the purposes of 

reconstruction. By means of the parallel rebinning, the measured values are resorted and 
reinterpolated as though they had been measured using a different radiation source (an 
expanded radiation source which is arranged on part of a helix and can emit beam fans that 
are in each case parallel to one another) and using a different detector (a flat, rectangular 

30 "virtual detector" including the axis of rotation 14). 

This is explained in more detail with reference to figure 3. Here, 17 designates 
the helix-like trajectory from which the radiation source irradiates the examination area. A 
fan-shaped beam cluster 43, the beams of which run in a plane that includes the axis of 
rotation 14, is emitted from the radiation source position S 0 . It is possible to consider that the 
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cone-shaped beam cluster emitted from the radiation source at position S 0 is composed of a 
large number of planar beam fans which are located in planes parallel to the axis of rotation 
14 and intersect at the radiation source position S 0 . figure 3 shows just one of these beam 
fans, namely the beam fan 43. 
5 Moreover, further beam fans 41, 42, and 44, 45 are shown in figure 3, said 

beam fans being parallel to the beam fan 43 and lying in planes that are parallel to one 
another and to the axis of rotation 14. The associated radiation source positions S_ 2 , S_ x and 
S l9 S 2 are taken up by the radiation source S respectively before and after it has reached the 
radiation source position S 0 . 

10 The beam fans 41 to 45 form a group and define a beam cluster 70 with a tent- 

like shape. A group of beam fans is called a projection. For each projection there is then 
defined a rectangular, virtual detector 160 which lies in a plane that includes the axis of 
rotation 14 and is oriented perpendicular to the parallel beam fans of a projection. The corner 
points of the virtual detector 1 60 are the penetration points of the beams which from the outer 

15 radiation source positions come into contact with the opposite helix section, through this 
plane. For the beam cluster 70 in figure 3, S_ 2 and S 2 are the outer radiation source 
positions. Detector elements that are arranged in a Cartesian manner are defined on the 
rectangular detector 160, that is to say rows and columns on which the measured values are 
reinterpolated. 

20 Subsequently, in step 107, the measured values allocated to the individual 

beams are multiplied by a weighting factor which corresponds to the cosine of the cone angle 
of the respective beam. The cone angle of a beam is the angle encompassed by this beam 
with a plane oriented perpendicular to the axis of rotation 14. If said angle is small, then the 
cosine of the angle is essentially equal to 1, so that step 107 may be omitted. 

25 In step 1 09, a one-dimensional filtering with a transmission factor that 

increases in a ramp-like manner with the spatial frequency is applied to the measured values. 
For this, use is made of respectively successive values in the direction perpendicular to the 
axis of rotation 14, that is to say along a row of the virtual detector 160. This filtering is 
carried out along each row of the virtual detector for all groups of beam fans. 

30 In other embodiments the parallel rebinning could be omitted. It is then known 

to modify the filtering since the detector unit is curved for example in the manner of an arc 
around the radiation source or around the axis of rotation. 
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In step 1 1 1 , the filtered measured values are then used for reconstruction of 
the absorption distribution in the examination area by means of a back-projection. The 
individual steps of the back-projection are shown in figure 4. 

In step 201 a voxel V(x) is defined within a predefmable field of view (FOV). 

5 Since the reconstruction is to take place with a low resolution, the amount of voxels may be 
for example 64 3 and the FOV may be 20 * 20 x 20 cm 3 . In step 203, a projection, that is to 
say a group of beam fans, is then selected which has not yet been used to reconstruct the 
voxel V(x) . If no beam of the projection runs centrally through the voxel F(x) , then the 
point at which a central beam would have come into contact with the detector surface is 

10 determined. The associated measured value is then calculated by interpolating the measured 
values of adjacent beams. The measured value which can be allocated to the beam of the 
projection that passes the voxel, or the corresponding measured value obtained by 
interpolation, is accumulated on the voxel V(x) in step 205. In step 207 a check is made as to 
whether all projections have been considered. If this is not the case, then the flowchart 

1 5 branches to step 203. Otherwise a check is made in step 209 as to whether all voxels V(x) in 
the FOV have been passed through. If this is not the case, then the method continues with 
step 201. If, on the other hand, all voxels V(x) in the FOV have been passed through, then 
the absorption in the entire FOV is determined and this reconstruction method is terminated. 

In order to determine within the examination area the spatial area in which the 

20 heart is located, in the following step 1 13 the heart is segmented in the reconstructed three- 
dimensional data record. 

A simple possibility for the segmentation is manual segmentation. A user, for 
example a physician, places markings on the surface of the heart in the three-dimensional 
data record. These markings are connected by lines and thus form a network that represents 

25 the surface of the heart. 

A further possibility for the segmentation consists in the use of a deformable 
model of the heart which is moved, rotated and scaled in the data record such that the 
correlation between the data record and the model of the heart is maximized. This known 
segmentation method is explained for example in "Deformable models in medical image 

30 analysis: A survey", Medical Image Analysis, 1(2): 91-108, 1996, to which reference is 
hereby made. 

Another known segmentation possibility is based on a region growing or 
region expansion process, in which a user predefines a so-called seed voxel in an object that 
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is to be segmented, that is to say in this case in the heart. Neighboring voxels of the seed 
voxel are then examined using an association criterion to determine whether they do or do 
not belong to the heart. This association criterion may be for example the fact of whether or 
not they are included within a range of values of the data values of the voxel. If a data value 
5 lies within the range of values then the corresponding voxel is assigned to the heart. In the 
next step, the neighboring voxels of the voxel which has been newly assigned to the heart are 
examined with respect to the association criterion and in some circumstances also assigned to 
the heart. This method is repeated until no more neighboring voxels can be assigned to the 
heart. 

1 0 The segmented heart shows the spatial area taken up by the heart in the 

examination area. 

This determination of the spatial area taken up by the periodically moving 
object in the examination area represents only one embodiment. Any method which makes it 
possible to determine this spatial area may be used according to the invention. For instance, 

1 5 the examination area could also be reconstructed using other known reconstruction 

techniques, such as the w-PI method which is described in "The «-PI Method for Helical 
Cone-Beam CT", R. Proska, Th. Kohler, M. Grass, J. Timmer, IEEE Transactions on 
Medical Imaging, Vol. 19, 848-863, September 2000. The reconstruction could also be 
carried out using the method from "A General Cone-Beam Reconstruction Algorithm", G. 

20 Wang, T.H. Lin, P.C. Cheng, D.M. Shinozaki, IEEE Transactions on Medical Imaging, Vol. 
12, 486-496, March 1993. The use of the reconstruction method known as ASSR (Advanced 
Single-Slice Rebinning) would also be possible, which although leading to poorer-quality 
images compared to the abovementioned reconstruction does require less computational 
complexity. This method is published for example in "Advanced Single-Slice Rebinning in 

25 Cone-Beam Spiral CT", M. KachelrieB, S. Schaller, W.A. Kalender, Medical Physics, Vol. 
27, No. 4, 754-772, 2000. A reconstruction method for a circular trajectory, along which the 
radiation source moves relative to the object, is shown for example in "3D Cone-Beam CT 
Reconstruction for Circular Trajectories", M. Grass, Th. Kohler, R. Proska, Physics in 
Medicine and Biology, Vol. 45, No. 2, 329-347, 2000. In reconstruction to determine the 

30 spatial area taken up by the moving object, it is important that the resolution of the 

reconstructed three-dimensional data record is selected such that the spatial contents can be 
determined for example by means of known segmentation methods. This condition is met, for 
example, when a volume of 20 * 20 * 20 cm 3 is represented by 64 3 or more voxels. Any 
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method capable of segmenting an object in a three-dimensional data record can be used 
according to the invention for the segmentation. 

Once the heart has been segmented in the three-dimensional data record, the 
geometric center of the heart can be determined in step 1 15 by simple geometric 
5 considerations. Only the geometric center of a known three-dimensional object is to be 

determined. The geometric center would be, for example, the center of gravity of this object 
if a spatially constant density of the object were to be assumed. 

In step 1 17, a point in time is allocated to each projection determined in step 
105. For this, it is determined which beam fan of a projection irradiates the geometric center 

10 of the heart. The point in time at which the measured values corresponding to this beam fan 
were acquired is allocated to the respective projection. 

If in other embodiments no parallel rebinning according to step 105 has been 
carried out during the determination of the spatial area taken up by the periodically moving : ,- 
object in the examination area, then said parallel rebinning must be carried out prior to step ■%£ 

15 117. tff 
For the reconstruction of the periodically moving object while taking into 
account the movement recorded during the measurement, repeating time ranges are selected 
in step 119, where the object has assumed at least a substantially identical object state in each 
time range. The at least one substantially identical object state thus occurs in each time range. 3- 

20 Time ranges in which the object states assumed by the object during these time ranges differ— 
from one another as little as possible are preferably selected. For example, two object states 
differ little from one another when the difference between two measured signals that 
characterize the respective object states, that is to say for example in the case of the heart the 
difference in the signals from the electrocardiograph, is small compared to the difference 

25 between the maximum and minimum value of the measured signal detected during a 

measurement. A difference between two measured signals is small, for example, if it is less 
than 1%, 2% or 5% of the difference between the maximum and minimum value of the 
measured signal detected during a measurement. In the reconstruction explained further 
below, use is only made of measured values which were acquired during these time ranges. 

30 In the case of the heart it is useful to select a time range of the diastolic phase of the 

movement of the heart, since in said phase the movement of the heart is considerably smaller 
than in the systolic phase. This is shown in figure 5. The period of the electrocardiogram H 
consists of a range Hi with relatively little movement and of a range H 2 with a lot of 
movement. In this example of embodiment the time range Hi is selected for the further 
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reconstruction since in this range the object states differ less from one another than in the 
range H 2 . 

In step 121, those projections which were parallel rebinned in step 105 and 
which were acquired during the time ranges Hi are determined. This is again shown in figure 
5. In step 117, points in time t 0 and ti were allocated to the projections P 0 and Pi. The points 
in time t 0 and ti lie in the time range Hi, so that the projections P 0 and Pi are used for the 
subsequent reconstruction. By contrast, projections to which the points in time t 2 and t 3 were 
allocated are not taken into account since the points in time t 2 and t 3 lie in the time range H 2 . 

If in other embodiments the measured values have not been multiplied 
according to step 107 by the cosine of the cone angle of the beam corresponding to the 
respective measured value during the determination of the spatial area taken up by the object 
in the examination area, then this can be carried out prior to the subsequent back-projection. 
The same applies to the filtering according to step 109. 

The measured values of the projections determined in step 121 are then used in 
step 123 for the reconstruction of the absorption distribution in the heart by means of a back- 
projection. The individual steps of this reconstruction are shown in figure 6. 

For this, in step 301a voxel V(x) is determined within a predefinable field of 
view (FOV). The amount of voxels may be for example 512 3 and the FOV may be 20 x 20 x 
20 cm 3 . In step 303, a projection, that is to say a group of beam fans, is then selected from the 
projections determined in step 121, which projection has not yet been used to reconstruct the 
voxel V(x) . If no beam of the projection runs centrally through the voxel V(x) , then the 
point at which a central beam would have come into contact with the detector surface is 
determined. The associated measured value is then calculated by interpolating the measured 
values of adjacent beams. The measured value which can be allocated to the beam of the 
projection that passes the voxel, or the corresponding measured value obtained by 
interpolation, is multiplied in step 305 by a weighting factor w p (x) . The subscript index in 
this case denotes the projection selected in step 303. 

The following consideration is used to determine the weighting factor w p (x) . 
During the acquisition, the voxel F(x) is irradiated for the first time by the cone-shaped 
beam cluster at a specific point in time. This point in time is referred to as the sunrise SR. 
The beam leaves the beam cluster again at a second point in time. This point in time is 
referred to as the sunset SS. A factor a p (x) is then allocated to the projection selected in 
step 303. This factor is smaller for projections whose point in time allocated in step 117 lies 
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relatively close to SR or SS than for projections whose corresponding point in time does not 
lie relatively close to SR or SS. Relatively close may mean for example that the point in time 
of the projection is located in a range that adjoins SR or SS, which corresponds to 5%, 10%, 
15% or 20% of the overall range between SR and SS. One exemplary profile of the factor 
5 a p (x) as a function of the position of the point in time of the projection between sunrise and 

sunset is shown in figure 7. In said figure, the time range between SR and SS is divided into 
three sections. In the first section Ti, which takes up 1 0% of the time range between SR and 
SS, the factor a p (x) increases from SR up to a value a c . In the following range T 2 , the 

factor is constantly equal to a c and in the third range T3, which also takes up 10% of the 

10 time range between SR and SS, the factor decreases again toward SS to zero. 

It is furthermore assumed that the movement of the heart is less pronounced at 
a point in time which lies in the center of the time range Hi, which comprises little 
movement, than at a point in time lying closer to the boundary with the time range H2, which 
comprises a lot of movement. Therefore, measured values of projections whose points in time 

15 lie in the center of the range Hi are more highly weighted than measured values of 

projections whose points in time lie closer to the boundary with the time range H2. In order to 
take this into account, a second factor P p (x) is assigned to the measured value. The 

subscript index p in this case denotes a projection and, since a point in time is allocated to 
each projection, also a point in time. Depending on the position of the point in time within 

20 the time range Hi, this factor may have a profile that ensures the following: the factor is 

greater for measured values of projections whose points in time lie in the center of the time 
range Hi than for measured values of projections whose points in time lie closer to the 
boundary with the time range H2. An exemplary profile is shown in figure 8. 

Since the same point in time has been allocated in step 1 17 to all measured 

25 values whose beams belong to a projection, all measured values of a projection have the same 
weighting factor w p (x) for a voxel V(x) , which weighting factor is defined by the following 

equation: 

2>/(x)A0O 

Here, (x)/?, (x) denotes a sum of all projections which have been determined in step 
30 121 and are not redundant. 
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In step 307, the weighted measured value is accumulated on the voxel V(x) . 
In step 309 a check is made as to whether all projections have been considered. If this is not 
the case, then the flowchart branches to step 303. Otherwise a check is made in step 311 as to 
whether all voxels V(x) in the FOV have been passed through. If this is not the case, then the 
5 method continues with step 301. If, on the other hand, all voxels F(x) in the FOV have been 
passed through, then the absorption in the entire FOV is determined and the reconstruction 
method is terminated (step 313). 

In other embodiments, step 123 and steps 301 to 3 13 may be replaced by other 
known reconstruction methods which generate a three-dimensional data record from the 
10 parallel projections determined in step 121 . 

Furthermore, the reconstruction in steps 301 to 3 13 may be restricted to 
projections which are not redundant. A projection is redundant with respect to a voxel V(x) 
if the beam of the projection which runs through this voxel V(x) is oriented in the opposite 
direction with respect to a beam of another projection that has already been used to 
15 reconstruct this voxel. 



WO 2004/089217 

CLAIMS: 



PCT/IB2004/001061 



1 . A computer tomography method comprising the steps: 

using a radiation source (S) to generate a cone-shaped beam cluster (4) which 

passes through an examination area (13) and a periodically moving object which is located in 

the examination area (13), 
5 - producing a relative movement between the radiation source (S) on the one 

hand and the object located in the examination area (1 3) on the other hand, where a 

trajectory, along which the radiation source moves relative to the object, runs on an 

imaginary cylindrical surface that surrounds the object, 

using a detector unit (1 6) to acquire measured values which depend on the 
10 intensity in the beam cluster (4) on the other side of the object, during the relative movement, 
recording the periodic movement of the object during the acquisition, 
reconstructing a spatial distribution of the absorption of the periodically 

moving object from the measured values with the aid of the recorded periodic movement of 

the object, comprising the steps: 
15 a) determining the spatial area taken up by the object in the examination area 

(13), 

b) subjecting the measured values to parallel rebinning in order to form a number 

of groups, where the beams corresponding to the measured values of each group form beam 
fans (41 . . .45) which lie in planes that are parallel to one another and to the axis of rotation, 
20 c) determining for each group a measured value whose beam irradiates the spatial 

area taken up by the object, and allocating to the respective group the point in time at which 
this measured value was acquired, 

d) determining those groups whose points in time, allocated in step c), lie within 

periodic, predefined time ranges (Hi), 
25 e) reconstructing the absorption distribution in the object from the measured 

values belonging to the groups determined in step d). 
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2. A computer tomography method as claimed in claim 1 , characterized in that 

the determination of the spatial area taken up by the object, in step a), comprises the 
following steps: 

reconstructing from the measured values a three-dimensional data record 
5 which contains the object, with a resolution which makes it possible to segment the object in 
the three-dimensional data record, 

segmenting the object in the three-dimensional data record, where the 
segmented object shows the spatial area taken up by the object in the examination area (13). 

10 3. A computer tomography method as claimed in claim 1, characterized in that in 

step c) the geometric center of the spatial area taken up by the object in the examination area 
is determined and for each group a measured value is determined whose beam runs through 
the geometric center, where the point in time at which this measured value was acquired is 
allocated to the respective group. 

15 

4. A computer tomography method as claimed in claim 1 , characterized in that 

the periodically moving object is a heart, where the periodic time ranges (Hi) are predefined 
with the aid of an electrocardiograph (8). 

20 5. A computer tomography method as claimed in claim 1, characterized in that 

the object moves less in the periodic, predefined time ranges (Hi) than in other time ranges 
(H 2 ). 

6. A computer tomography method as claimed in claim 1, characterized in that 
25 the reconstruction is carried out with the aid of a filtered back-projection. 

7. A computer tomography method as claimed in claim 1, characterized in that 
the relative movement between the radiation source (S) on the one hand and the object 
located in the examination area (13) on the other hand comprises a rotation about an axis of 

30 rotation (14) and runs in a circular or helix-like manner. 

8. A computer tomography scanner, in particular for carrying out the method as 
claimed in claim 1, comprising: 
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a radiation source (S) for generating a cone-shaped beam cluster (4) which 
passes through an examination area (13) and a periodically moving object which is located 
therein, 

a drive arrangement (2, 5) for rotating the object located in the examination 
5 area and the radiation source (S) relative to one another about an axis of rotation (14) and 
moving them relative to one another parallel to the axis of rotation (14), 

a detector unit (16) for acquiring measured values, said detector unit being 
coupled to the radiation source (S) 5 

a movement recording device (8), in particular an electrocardiograph, for 
10 recording the periodic movement of the object during the acquisition, 

at least one reconstruction and image processing computer (10) for 
reconstructing the spatial distribution of the absorption within the examination area (13) from 
the measured values acquired by the detector unit, with the aid of the periodic movement of 
the object recorded by the movement recording device (8), 
1 5 a control unit (7) for controlling the radiation source (S), the drive 

arrangement (2, 5), the detector unit (16), the movement recording device (8) and the at least 
one reconstruction and image processing computer (10) in accordance with the following 
steps: 

using a radiation source (S) to generate a cone-shaped beam cluster (4) which 
20 passes through an examination area (13) and a periodically moving object which is located in - 
the examination area (13), 

producing a relative movement between the radiation source (S) on the one 
hand and the object located in the examination area (13) on the other hand, where a 
trajectory, along which the radiation source moves relative to the object, runs on an 
25 imaginary cylindrical surface that surrounds the object, 

using a detector unit (1 6) to acquire measured values which depend on the 
intensity in the beam cluster (4) on the other side of the object, during the relative movement, 
recording the periodic movement of the object during the acquisition, 
reconstructing a spatial distribution of the absorption of the periodically 
30 moving object from the measured values with the aid of the recorded periodic movement of 
the object, comprising the steps: 

a) determining the spatial area taken up by the object in the examination area 

(13), 
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b) subjecting the measured values to parallel rebinning in order to form a number 
of groups, where the beams corresponding to the measured values of each group form beam 
fans (41 . . .45) which lie in planes that are parallel to one another and to the axis of rotation, 

c) determining for each group a measured value whose beam irradiates the spatial 
5 area taken up by the object, and allocating to the respective group the point in time at which 

this measured value was acquired, 

d) determining those groups whose points in time, allocated in step c), lie within 
periodic, predefined time ranges (Hi), 

e) reconstructing the absorption distribution in the object from the measured 
10 values belonging to the groups determined in step d). 

9. A computer program for a control unit (7) for controlling a radiation source 

(S), a drive arrangement (2, 5), a detector unit (16), a movement recording device (8) and at 
least one reconstruction and image processing computer (10) of a computer tomography 
15 scanner, in particular for carrying out the method as claimed in claim 1, as follows: 

using a radiation source (S) to generate a cone-shaped beam cluster (4) which 
passes through an examination area (13) and a periodically moving object which is located in 
the examination area (13), 

producing a relative movement between the radiation source (S) on the one 
20 hand and the object located in the examination area (13) on the other hand, where a 
trajectory, along which the radiation source moves relative to the object, runs on an 
imaginary cylindrical surface that surrounds the object, 

using a detector unit (16) to acquire measured values which depend on the 
intensity in the beam cluster (4) on the other side of the object, during the relative movement, 
25 - recording the periodic movement of the object during the acquisition, 

reconstructing a spatial distribution of the absorption of the periodically 
moving object from the measured values with the aid of the recorded periodic movement of 
the object, comprising the steps: 

a) determining the spatial area taken up by the object in the examination area 
30 (13), 

b) subjecting the measured values to parallel rebinning in order to form a number 
of groups, where the beams corresponding to the measured values of each group form beam 
fans (41 . . .45) which lie in planes that are parallel to one another and to the axis of rotation, 
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c) determining for each group a measured value whose beam irradiates the spatial 
area taken up by the object, and allocating to the respective group the point in time at which 
this measured value was acquired, 

d) determining those groups whose points in time, allocated in step c), lie within 
periodic, predefined time ranges (Hi), 

e) reconstructing the absorption distribution in the object from the measured 
values belonging to the groups determined in step d). 
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